Field-induced Berezinskii-Kosterlitz-Thouless transition and string-density plateau in 
the anisotropic triangular antiferromagnetic Ising model 
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The field- induced Berezinskii-Kosterlitz-Thouless (BKT) transition in the ground state of the 
triangular antiferromagnetic Ising model is studied by the level-spectroscopy method. We analyze 
dimensions of operators around the BKT line, and estimate the BKT point He — 0.5229 ± 0.001, 
which is followed by the level-consistency check to demonstrate the accuracy of our estimate. Further 
we investigate the anisotropic case to clarify the stability of the field-induced string-density plateau 
against an incommensurate liquid state by the density-matrix renormalization-group method. 

PACS numbers: 64.60.-i, 05.50.+q, 05.70.Jk 



It was exactly proven that the nearest-neighbor (NN) 
antiferromagnetic Ising model on the triangular lattice 
shows no phase transition at finite temperature, and pos- 
sesses a ground-state ensemble carrying the finite residual 
entropy per spin, S ~ 0.3231fcB This circumstance 
stems from the frustration effect that not all spins at the 
corners of each elementary triangle can be energetically 
satisfied, and then brings about the power-law decays 
of the correlation functions of physical quantities in the 
ground state ■ There is a long history of research on 
various types ofperturbation effects on this ground-state 
degeneracy 0, 0, H, IE S Uli where, as we will see in 
the following, the exact mapping to the so-called trian- 
gular Ising solid-on-sohd (TISOS) model 0,113 followed 
by a coarse graining provides an effective field theory 
to describe the low-energy physics addition, it 

should be remarked that the ground-state spin configura- 
tions can be classified according to the number of strings 
(see below) , which provides an intuitive connection to 
one-dimensional (ID) quantum systems with global U(I) 
symmetry under the path-integral representation. 

In this paper, we treat an anisotropic triangular an- 
tiferromagnetic Ising model (TAFIM) under a magnetic 
field; its reduced Hamiltonian H = /JHtafim is given as 




(1) 



The binary variable aj ~ 0, 1 is on the jth site of the 
triangular lattice A which consists of interpenetrating 
three sublattices A; {I = 0,1,2), and the first (second) 
sum runs over all NN pairs (sites). The AF coupling Kjk 
takes two values Ki -f /i or Ki depending on whether 
the bond (j, k) lies in the xi direction or not (see Fig. 

Here we define a quantity Q = ^jid2 with 
N-j, = 1 — Sa ■ o- j_, for all 79, and further restrict 
ourselves to the zero temperature case Ki — > 00. Then Q 
is independent of j2 and the Boltzmann weight per row 
is given by e'*'^. Therefore, the anisotropy parameter /i 
plays a role of the chemical potential to control the num- 
ber of strings (an example of the string representation is 
given in Fig.^ 0. Our main goal is to clarify the phase 



FIG. 1: An example spin configuration with the \/3 x \/3 
structure. The jih site is specified by two integers (^1,^2) as 
labeled. The long (short) side of the rectangle frame is in the 
xi (2:2) direction. The spins on two (one) of three sublattices 
Ao and Ai (A2) are parallel ® (antiparallel (g)) to the field 
direction. Four gray stripes in the X2 direction give the string 
representation of the spin configuration. 



diagram of the model ^ in its ground state Ki — > 00. 
For this, we shall use the level-spectroscopy (LS) method 
[T^ to treat the Berezinskii-Kosterlitz-Thouless (BKT) 
transition induced by H in the isotropic case = 0. 
On the other hand, for the anisotropic case ^ 7^ 0, the 
Pokrovski-Talapov (FT) transition [13 between a com- 
mensurate ordered phase and an incommensurate liq- 
uid phase is expected. We directly calculate the ^ de- 
pendence of the number of strings by the use of the 
density-matrix renormalization-group (DMRG) method 
|l4j . Then, we provide a reliable phase diagram. 

As we will see in the following, the string degrees of 
freedom in which the frustration effects arc encoded play 
a central role for both understanding the phase diagram 
and relating it to the magnetization plateaus observed 
in the ID frustrated quantum spin systems |l5l |. This 
comes from the fact that the quantity Q corresponds to 
the uniform magnetization, and thus fi can be regarded 
as the magnetic field in the quantum spin systems. Fur- 
ther, quite recently, the string representation has been 
also employed in investigations of the 2D frustrated quan- 
tum magnets [T^ . so the precise analysis on the present 
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TABLE I: Estimates of He- We label the estimate by 
the spin-wave excitation (the vortex excitation) as KTl 
(KT2). We also give the results by the phenomenological 
renormalization-group (PRG) and the level-spectroscopy (LS) 
methods. 



Reference 


Method 


Extrapolation 




Error 


m 


KTl 


A rough estimate 


0.6 




m 


KTl 


Iterated fit 


0.532 


±0.02 




KT2 


Iterated fit 


0.52 


±0.04 


m 


PRG 


Finite-size scaling 


0.422 


±0.014 


This work 


LS 


Polynomial fit 


0.5229 


±0.001 



fundamental model may also offer a hint for a further 
understanding of quantum systems. 

Now we shall start with the isotropic case /i = 0. 
For H = 0, from the exact asymptotic behavior of the 
spin-spin correlation function the scaling dimension of 
the staggered magnetization (S) is given hy xs = j Q, 
while the dimension of the uniform magnetization (s) is 
hy Xs = J 01- Thus the magnetic field is irrelevant, and 
the critical region continues up to a certain value H^. For 
H > He, the criticality of the ensemble disappears and 
the threefold-degenerate ground state with the -s/S x V3 
structure of the sublattice is realized, where a majority 
spin is in the field direction (see Fig.^. The transition 
at He is the BKT type, and is described, in the scaling 
limit, by the 2D sine-Gordon Lagrangian density ^3 



2ttK 



■ cos 3V2(f), K 



(2) 



where y cx H and the continuous field in the 2D Eu- 
clidean space {xi,X2) proportional to the height vari- 
able of the TISOS model satisfies V2(l> + 2tt = V^cj) [13 ■ 
Since the lattice model is not exactly solvable, there have 
been several attempts to numerically estimate He'- Blote 
and Nightingale investigated this problem in detail by 
the transfer-matrix method 0, Q ■ Actually, they eval- 
uated finite-size estimates Hc{L) by numerically solving 
the equation for the scaled gap, i.e., the so-called KT 
criterion xs{H,L) = | (this is referred to as KTl in 
Table P) . Then, in order to accelerate the slow conver- 
gence of Hc{L), the iterated fits with taking account of 
the logarithmic correction were performed (for a more re- 
cent estimation, see 0). On the other hand, de Queiroz 
et al. treated the same model by the phenomenological 
renormalization-group (PRG) method [H; they exhib- 
ited a much smaller value inconsistent with the previous 
estimations (see Ref. 0). However, it is often pointed out 
that PRG calculations fail to estimate the BKT points 
|l9j . so their result may suffer from an inadequacy of 
the method. Consequently, to accurately determine H^, 
there still remains some difficulty. 

In the studies of ID quantum systems, however, the 
LS method provides an efficient way to treat the BKT 
transitions [T^l- This is also true for the investigations 
of 2D classical spin systems |23|. Let us consider the 
system on A with M (— s- oo) rows in the X2 direction of 



L (a multiple of 3) sites in the xi direction wrapped on 
the cylinder and define the transfer matrix connecting 
the next-nearest-neighbor rows. Since, in our discussion, 
the number of the strings Q (or its density p = Q/L) 
is the most important conserved quantity in the trans- 
fer, we explicitly specify a block of the matrix as Tq(L) 
and denote its eigenvalues as Xp^Q{L) or their logarithms 
as Ep^q{L) — — i In |Ap,Q(i)| {p specifies a level). In 
the isotropic case, the smallest one corresponding to the 
ground state is in the block Qo = 2L/3 nil; we shall 
denote it and the excitation gaps from it as -Bg.Qo {L) 



and AEp^Q{L) 



Ep.Q W ~ -^g,Qo (^) ' respectively. Then 



the conformal invariance provides direct expressions for 
the central charge c and a scaling dimension Xp,Q in 
the critical system as Eg^Qg{L) ~ Lf — ttc/QLC, and 
AEp^Q{L) ~ 2ttXp^q/LC. Here C 2/V3) and / are 
the geometric factor for A and a free energy per site, 
respectively jHEl. 

Blote and Nightingale precisely checked various scal- 
ing dimensions based on the Coulomb-gas scenario Q, 
whereas Nomura pointed out the importance of loga- 
rithmic corrections in the renormalized scaling dimen- 
sions x{l) = AE(L)/(2tt/LC) to determine the BKT 
point {I = hiL) [ij. In the present effective theory ||2Jl, 
there are two marginal operators on the BKT point, i.e., 
M = (l/K) (V0) and s = V^cosSV^cf) which hybridize 
along the RG flow and result in two orthogonalized ones, 
i.e., the "A^-like" and the "cos-like" operators Writ- 
ing the former and the latter as Oq oc (A^ + s/\/2) and 
Oi (X (— / + s) , their renormalized scaling dimen- 
sions can be calculated near the multicritical point; the 
results up to the first-order perturbations are given as fol- 
lows: a;o(0-2-yo(l + |i) and xi{l) ~ 2 + yo {2 + ^t) , 
where (j/oj^i) = (9/\/2 — 2,y) and the small deviation 
from the BKT point t = yi/yo — 1 |l2l |. Another impor- 
tant operator is a relevant one, i.e., the staggered magne- 
tization S = exp(±i'\/20) whose dimension is expressed 
as xs{l) — I (l ± ■^yo) in the same region. Consequently, 
the level-crossing condition. 



xo(/) = 4-9xs(0, 



(3) 



provides a finite-size estimate of the BKT point. Since 
these operators are described by 0, we can calculate 
the renormalized scaling dimensions from excitation gaps 
found in the Qg block. Further, the symmetry proper- 
ties such as the translation of one lattice spacing (or a 
cyclic permutation among sublattices A^) T, the space 
inversion V, and the spin reversal 5 at _ff = are also 
important for the precise specification of relevant exci- 
tations levels. These symmetry operations can be inter- 
preted in the field language as T: ^/2(f>+2TT/3, V: 
V2(t> ^ -V20, and S: V2<^ ^ ^/20 + tt |1 US . For 
instance, s is invariant for T and V, but is odd for S as 
expected, while S is in the = 27r/3 block and is odd 
for S. Therefore, the corresponding excitations to Oq 
and Oi can be found in the subspace of the wavevector 
k = and the even parity for the space inversion. We cal- 
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H=0, Exact: arccos(1/2' 
H=2, DMRG L=60 
H=4, DMRG L=60 



FIG. 2: The spectroscopy of levels observed in the L = 24 
sites system. The crossing of a;o(/) (diamonds) and 4 — 9a;g(Z) 
(triangles) gives the finite-size estimate of the transition point 
Hc{L). Inset shows the extrapolation of finite-size estimates 
to the thermodynamic limit L — > oo. 

TABLE II: The L dependences of dimensions xo{l) and xi{l) 
at He- The average x^^{l) (see text) is extrapolated to L ^ cxj 
using the least-squares fitting of the polynomial in Ijl? . 

L 18 21 24 27 30 oo 

xo{l) 1.88528 1.88110 1.87913 1.87835 1.87823 ~ 

xiil) 2.40464 2.37543 2.35551 2.34096 2.32981 — 

x^Sjl) 2.05840 2.04587 2.03792 2.03256 2.02875 2.01292 



culate the excitation levels by utilizing these symmetry 
operations. 

We perform the exact-diagonalization (ED) calcula- 
tions of Tqq(L) for systems up to L = 30. We show 
the level-crossing data in Fig. |21 and the extrapolation 
of the finite-size estimates Hc{L) to the thermodynamic 
limit L — > oo using the least-s quar es fitting of the poly- 
nomial in (the inset) Then we find that 
while our result is consistent with previous estimations, 
it is much more accurate owing to the fast convergence of 
finite-size estimates (see Tabled). Next we shall check a 
universal relation among excitation levels. For instance, 
the following relation is to be satisfied at the BKT point: 
[2a;o(/)+a;i(/)]/3 ~ 2 [13. In Table|Ill we give the scaling 
dimensions at i?c, where the left-hand side of the relation 
is denoted as Xa.-v{l)- Although xo{l) and xi{l) consider- 
ably deviate from the value for the free-boson case 2 due 
to the logarithmic corrections, their main parts cancel 
each other, so the average takes a value close to 2. These 
data provide the check of the accuracy of our estimate 
and the evidences to ensure that numerically studied lev- 
els have the above-mentioned theoretical interpretations. 

In the rest part, we shall discuss the anisotropic case 
/i 7^ 0. As we have seen in the above, the magnetic field 
H favors the long-range ordered commensurate (C) phase 
through the potential cos3-\/2^. On the other hand, ^ 
which newly introduces a local density term di(j) to the 
effective theory tends to stabilize an incommensurate 
(IC) liquid phase 0, [2^. Therefore, as can be found 



!-*'-1)/7t 

0^2 4 

FIG. 3: The ^-p curves: Triangles (diamonds) with dot- 
ted fine show the DMRG data at = 2 (i? = 4), and the 
solid line exhibits the exact result aX H = Q. The flat re- 
gion with p = I corresponds to the string-density plateau 
[^-{H), yL+{H)\. /io(-ff) is the threshold below which the 
string is absent. 



in the literature, the PT-type C-IC transition may oc- 
cur 0, H^. For _ff = 0, we know the exact ^x-p curve 
p(/i) = arccos (l/2e^'^ — l) /tt which is given in Fig. |3| 
01 . On the other hand, for H ^ 0, we employ numeri- 
cal methods and estimate the curve from the finite-size 
system data as ^ ~ [Eg^Q+2iL) — Eg^Q{L)] /2 @. While, 
like the exact solid curve, p is a smooth function of fi 
showing the compressive liquid state for H < He, there 
is the string-density plateau [/i_(77),/i_|_(if)] with p = | 
for H > He (see the DMRG data denoted by marks). 
In this plot, the top and the bottom of each step cor- 
respond to {Q + 2)/L and Q/L respectively, so we can 
estimate the C-IC phase boundary lines from the edges 
of the plateau. 

In Fig. 21 we provide our phase diagram. Here it is 
noted that the threshold fJ,o{H) below which the doubly 
degenerate vacuum of strings with p = is realized (see 
Fig El) is exactly given by p,a{H) = - ln[2 cosh(i7/4)] 
0, l27l |. So, we also draw it in the figure. The cross 
on the isotropic line /i = shows the BKT point ob- 
tained by the LS method. For large H, we can use the 
ED data fx±{H,L) ^ ± [Eg^Q,±2{L) - Eg^Q,{L)] /2 and 
the extrapolation formula ii± {H,L) ~ /.t± (iJ ) -I- const/L^ 
(see open circles) |2^. For H < 4, assuming the square- 
root behavior around the plateau, we estimate fi±{H) 
from the p-fi curve obtained by DMRG. From this plot, 
we find that two PT-transition lines )jl±{H) seem to be 
terminated at the BKT point (//, H) = (0, He) and that 
the plateau region becomes wider with the increase of 
the magnetic field. For i7 ~ iJc, it is still difficult to de- 
termine the narrow plateau region corresponding to the 
exponentially small energy gap even by the use of the 
DMRG method. However, by combining the LS result 
and the DMRG data, we can obtain the reliable phase 
diagram of our model. 

Lastly, we shall discuss some related topics. The mag- 
netization process observed in the ground state of the 
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FIG. 4: The ground-state phase diagram. The vacuum of 
strings is stabilized on the left of the solid curve iJ,o{H). The 
field-induced ordered phase corresponds to the string-density 
plateau state with p = |. The cross on the fj, = line shows 
the BKT point obtained by the LS method. The filled (open) 
circles show estimates ^±{H) by the DMRG (ED) method. 
Dotted curves give a guide to the eyes. 

S = anisotropic AF chain with a strong frustration 
exhibits the plateau at i of the saturation magnetiza- 
tion |l5j . This plateau state exhibits spontaneous break- 
ing of the translational symmetry down to the period 
n = 3, so it is threefold degenerate (i.e., TTiTTi ■■■): 
and thus satisfies the necessary condition for the magne- 
tization plateau |29l |. i.e., n{S — m) = integer with the 



average magnetization per site m = g . According to the 
bosonization treatment |3Cll | , this magnetization-plateau- 
formation transition is described by the sine-Gordon field 
theory, which is identical to the present case. And more 
generally, since Q is the conserved quantity in our sys- 
tem, the above necessary condition can be translated 
to the density plateau condition for the string systems 
as np = integer, where p is the possible density at the 
plateau and n is the periodicity of the plateau state. For 
its derivation, U = exp[— z ^^^^-^ (27rji/L) A'^j^ j-^] plays 
a role of the twist operator. While the present system 
only possesses the p = | plateau state with the spon- 
taneous symmetry breaking of T down to the period 
n = 3, Noh and Kim intensively investigated the in- 
teracting string/domain- wall systems based on TAFIM 
with spatially anisotropic further-neighbor couplings |2(il | . 
They used the Bethe-ansatz method to diagonalize the 
transfer matrix, and found the p = ^ plateau phase. 
Furthermore the period n = 2 of the plateau state has 
been suggested in that case. Therefore, we think that 
the above necessary condition for the plateau can pro- 
vide a important viewpoint for the understandings of the 
interacting strings embedded in TAFIM with various ex- 
tensions. 
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